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I. INTRODUCTION 

Compared with mechanisms for the growth of water waves, those 
for wave decay have received little attention. In wave forecasting, 
for example, there is far more agreement on empirical wave growth 
formulas than on decay formulas. In the case of decay due to ambi- 
ent ("external") underwater turbulence, the relative neglect is 
probably due to the lack of both ocean turbulence data and a simple 
theory. That of Phillips (1959) is based on the turbulence spec- 
trum near a free surface, which has not been studied in detail. 
Earlier, more speculative work (Groen, 1954; Bowden, 1950) is also 
closely tied to unknown characteristics of the near-surface turbulence 
However, wave decay is quite important, and the results of the 
simple but novel laboratory experiment described below provide some 
insight into this problem. These results are admittedly fragmentary: 
the experiment was preliminary in nature, and the work had to be 
abandoned unexpectedly. Nevertheless, it seems worthwhile to pre- 
sent our data and tentative conclusions to guide others in similar 
work . 

II. APPARATUS 



A paddle (or fan) generated surface waves which propagated 
down a long tank, over a zone of underwater turbulence, to a non- 
reflecting beach (Fig. 1). The turbulence was produced by a ver- 
tically oscillating underwater grid, which spanned the width of the 
tank, and consisted of l.-cm wide metal strips woven about 1 . -cm 
apart so that the energy containing eddies had a scale of 1. cm 
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(see Paquin, 1968, for details) . The oscillation frequency was 
approximately 1. Hz; the amplitude was 1. cm. The eddies had 
little direct effect on the water surface. The waves produced by 
turbulence alone were less than a millimeter in height, whereas 
the paddle- or wind-generated wave heights were on the order of 
centimeters upstream of (i.e., before entering) the turbulent 
region. No wave breaking was observed. 

Surface waves were shielded from mean currents caused by the 
grid drive mechanism by a false bottom (Fig. 1). Any such current 
remaining was completely masked by the velocity fluctuations, so 
that the effect of wave-current interactions was probably quite 
small. Wave frequencies ranged from 1 to 10 Hz. Frequencies above 
2 Hz corresponded to deep-water waves (Stoker, 1957), so that the 
influence of both bottom and grid was negligible. For lower fre- 
quencies the false bottom was removed (introducing weak currents 
not otherwise present) but the waves did feel the grid itself. 

Faraday waves (McGoldrick, 1968) were suppressed by thin metal 
strips near the paddle. 

Resistance wave gauges were mounted upstream and downstream 
of the turbulent zone, and were calibrated both statically and 
dynamically (using progressive waves and stroboscopic light). The 
response was flat from DC to 4 Hz and 10% down at 26 Hz. A con- 
stant-temperature anemometer (DISA Model 55A01) was used with a 
quartz-protected hot-film probe (DISA 55A81) to measure the turbu- 
lent motions. All data were FM recorded with a Precision Instruments 
6200 Tape Recorder and analyzed, after playback at 10 times or 100 
times the recording speed, with a General Radio Type 1900 Wave Analyzer. 
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MEASUREMENTS 



All measurements 


were 


made during 


steady condi 


tions with clean 


water. The decay 


due 


to viscosity 


and possible 


surface 


f i 1ms 


determined first, 


then 


( immediately 


afterward ) 


that due 


to these 



causes plus turbulence. 

A typical paddle-generated wave spectrum (with or without turbu- 
lence) is shown in Fig. 2. In order to investigate the widths of 
the harmonics, the recorder playback speed and the analyzer bandwidth 
were varied to produce an effective bandwidth that ranged from 1.0 Hz 
to 0.1 Hz. Since the observed energy level did not depend on the 
bandwidth, the wave energy was assumed to be located only at- the 
paddle frequency and its harmonics, and results were interpreted as 
line spectra. Figure 3 shows typical line spectra in the presence 
of turbulence . 

The turbulence velocity spectrum obtained with the anemometer 
placed in the center of the turbulent region at a depth of 2. cm, is 
shown in Fig. 4. Similar spectra were found at depths of 1. and 3. 
cm at various positions in the region. The lack of a large super- 
imposed unidirectional flow introduced serious problems associated 
with variations of probe sensitivity with flow direction. Calibra- 
tion showed the minimum sensitivity to be about 25% of the maximum. 
(The hot-film probe was always set to be most sensitive to flow in 
the wave propagation direction.) Thus, the turbulence spectrum 
should be interpreted in a qualitative sense, showing only that the 
spectral peak is below the wave frequencies, and giving some 
indication of the turbulence intensity. 
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A corresponding velocity spectrum due to paddle waves 
alone is also shown in Fig. 4. Despite the sensitivity prob- 
lem, the peaks associated with the paddle waves dominate/ and 
the energy associated with the fundamental wave is within a 
factor of two of that predicted by linear theory. 

Within the interacting region the turbulence energy is 
comparable to the wave energy. This is convenient for the 
experiment since large (hence easily measured) wave decay 
results. On the other hand/ it probably does not reflect con- 
ditions in a natural environment/ where wave motions are 
expected to dominate (Shonting, 1968) . 

IV. DISCUSSION 
A. Paddle-generated Waves 

According to laminar viscous theory for surface waves 
(Van Dorn/ 1966)/ wave decay with distance is given by the 
linear expression 



da/dx = 

which gives 

a ( x ) = a ( x q ) expt-g^ (x - x q )}. (1) 

Here, a is wave amplitude, x is propagation distance, and (3 , 

& 

the damping coefficient, varies with frequency f, viscosity, and 
the character of any surface film present. This is a first-order 
result in ka (k is wave number), valid for ka<<l. The measured 
ka values ranged from 0.04 to about 1; most were near 0.5. 
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Nevertheless, we use (1) to determine 8^. In no case was the 
observed laminar decay greater than could be explained by the 
presence of a surface film. 

Following Phillips (1959), two turbulent mechanisms can 

lead to wave decay. Scattering (due to turbulent convective 

distortion of wave fronts) is a first-order process in ka , 

results in the generation of a broad spectrum of scattered 

waves, and leads to an exponential decay law analogous to (1) 

with 8^ replaced by 8 fc - Phillips' results can be rewritten to 

show that two-dimensional waves passing through a horizontally 

homogeneous turbulent region experience scattering losses 
8/3 

giving 8 f 

Dissipation (whereby mean wave motions strain turbulent 
vortex filaments, thus increasing turbulent energy at the 
expense of wave energy) is a second-order process in ka , which 
is effective when the scale of the turbulence is much less than 
that of the waves. The decay law associated with this mechanism 
is not known, but is probably nonlinear in wave amplitude (in 
the differential form). 

It is likely that dissipation dominated in this experiment 
for the following reasons: (1) the mean square surface height 

downstream of the turbulent zone was typically less than 10% of 
that upstream. (2) Any decay due to scattering should be accom- 
panied by a broadening of the downstream wave spectrum; we could 
not detect any broadening with our filters. (This supported the 
observation that the downstream waves continued to appear long- 
crested.) (3) Finally, scattering should cause random fluctuations 



6 



in the contribution of any fixed frequency band to the rms sur- 
face height, as read on the wave analyzer voltmeter. (Note that 
the tank walls were smooth, allowing good reflection, so that 
wall effects on scattering should be slight.) Fluctuations of 
about 10 % of the mean value were observed downstream; they were 
undetectable upstream. Therefore, scattering accounted for a 
very small fraction of the total energy lost by the incident wave. 

B. Wind-generated Waves 

Wind waves generated in the laboratory obey an exponential 
growth law (Hidy and Plate, 1966) 

da/dx = g a . 

w 

Since both the growth due to wind and the laminar decay follow 

an exponential law, they can be combined to give an effective 

3 in ( 1 ) (which may no longer be positive). 

& 

The energy is now broad band. Measurements and calculations 
were done using energy in 1.- Hz bandwidths . The results are 
plotted against the center frequency of the filter band. 

IV. INTERPRETATION OF RESULTS 

The decay measurements are interpreted in two ways. 

A. Linear Decay Hypothesis 

Figure 5 shows the effect of assuming a linear turbulent 
decay process (such as would occur with scattering), which is 
independent of laminar decay: 
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da/dx = - ( 6^ + 6 fc ) a 
or 

a (x ) = a ( x ) exp {- ( 8 + 6 ) (x - x ) }. (2) 

O X, t o 

The observed decay with no turbulence gives ft (eqn. 1); the cor- 

A / 

responding decay with turbulence then gives 

The fundamentals travel as free waves. The second harmonics 
are not free waves, but bound waves associated with a periodic 
nons inusoidal profile (Lamb, Art. 250). We are interested in the 
decay of free waves. A reasonable correction is to plot the 
second- harmonic turbulent damping coefficients against the fre- 
quency of free waves of the same wave number. In the same spirit 
this correction is based on linear theory. (Aside from some 
third-order work by Wilton, 1915, there are no nonlinear disper- 
sion results for capillary-gravity waves.) A free wave obeys the 
dispersion law 



4 7 T 2f2 _ gk + Tk 3 /p (3) 

where g is gravity, p density, and T surface tension. Since the 
second-harmonic wave numbers and frequencies are double those for 
the first harmonic, the actual dispersion law for the second 
harmonics is 



7T 2 f 2 2 = gk 2 /2 + Tk 2 3 /8p 



Then the corrected second-harmonic frequency is 



(f 2 ) 



1 + Tk 2 2 /pgf /2 



cor r 



2 + Tk z /2Pg 
2 J 



* 2 - 



All second-harmonic damping coefficients are plotted at their 
corrected frequencies. 
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The paddle-wave results can be fitted to the power-law 



relation 




where is 3.2 x 10 sec cm (first harmonics), and c^ is 

- 3 2 -1 

1.8 x 10 sec cm (second harmonics). Based on the reasoning 
given above, we view the reasonable agreement with Phillips* 
8/3-power prediction as accidental. Also, Phillips* result rests 
on a model which is quite unlike the experiment: surface tension 

is neglected; the incident wave (already assumed small: ka<<l) 

is taken to vary little over the time of the decay calculation; 
the turbulence spectrum is fixed on dimensional reasoning based 
on the existence of an inertial subrange. All of these assump- 
tions are violated here. 

The discrepancy between values of c for the first two har- 
monics suggests either an energy transfer between the harmonics 
or, since the first harmonics are stronger than the second, that 
there is an amplitude dependence hidden in the damping coeffi- 
cient McGoldrick (1965) has shown that long-crested capillary- 

gravity waves do experience a resonant interaction between the 
first and second harmonics near the fundamental wave number 
K = (g/2T) (2.6 cm ^ for pure water at 20°C). For measured 
values of T, however, K was far above the experimental wave num- 
bers, so that resonance was highly unlikely. 
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B. Quadratic Decay Hypothesis 



The alternative explanation for the discrepancy between the 
first- and second-harmonic damping coefficients is explored by 
assuming (somewhat arbitrarily) a quadratic decay law, such as 
might be expected with turbulent dissipation, adding to the linear 
laminar decay: 



or 



da/dx = + Y fc a 2 ) 

a (x ) exp { - B . (x - x )} 
o £ o 

Y a(x ) 

1 + jj — [1 - exp {- 6^(x - x q ) }] 



(4) 



(5) 



The observed turbulent decay and eqn. 5 give y . The results 
are shown in Fig. 6. Although there is still some scatter, a 
more consistent picture has emerged. Aside from two questionable 
experimental points (common to both hypotheses) , all points can 
be fit to the power-law relation 

Y t = Cf 5 (6) 

- 5 5 - 2 

where C is 1.7 x 10 sec cm 

On the basis of this evidence, we suggest that steep capil- 
lary-gravity waves in the presence of underwater turbulence obey 
a dissipative decay law which is close to quadratic (i.e., 
eqn. 4) when 

(i) the length scale of the turbulence is much less than 
the wave length 
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(ii) the energy per volume of the turbulence is the same order 
as the maximum energy per volume of the waves. 
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